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The formation dynamics is studied for a singular profile of a surface of an ideal 
conducting fluid in an electric field. Self-similar solutions of electrohydrodynamic 
equations describing the fundamental process of formation of surface conic cusps with 
angles close to the Taylor cone angle 98.6° are obtained. The behavior of physical 
quantities (field strength, fluid velocity, surface curvature) near the singularity is 
established. 

It is known that a flat boundary of liquid metal becomes unstable in a strong elec- 

tric field. The development of instability results in conic cusp singularities, from which the 



strengthened field initiates emission processes 



The description of these processes 



is a key problem of the electrohydrodynamics of conducting fluids with free surfaces; inter- 
est in this problem is largely caused by the practical use of liquid-metal sources of charged 
particles. The progress in this field is associated with Taylor's work where it was demon- 
strated that the surface electrostatic pressure Pe for a cone with angle 98.6° depends on 
the distance from its axis as and, hence, can be counterbalanced by the surface pressure 
Ps ~ Since the force balance is violated at the cone apex, Taylor's solution cannot 
be treated as the exact solution of the problem of equilibrium configuration of a charged 
surface of conducting fluid and only represents the possible asymptotic form at r — cxd. 
At the same time, it turned out that Taylor's solution nicely describes the experimentally 
observed surface shape before the instant of singularity formation. It was pointed out in 

0; 0; ^ that thc auglc of incipient conic formations is close to the Taylor cone angle. 

What is the reason for such a coincidence? One may assume that the mechanism for the 
formation of conic cusps with an angle of 98.6° during a finite time is not directly associated 
with the static Taylor model. A high reproducibility of experimental results and a weak 
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dependence of fluid behavior at the final instability stages on the geometry of the system 
suggest that the behavior of fluid near the singularity has a self-similar character. 

Let us check the validity of this hypothesis. Consider the potential motion of an ideal 
fluid occupying the region bounded by free surface z = ri{x,y,t). We will assume that the 
vector of an external electric field is directed along the z axis and equals E. The velocity 
potential $ of fluid and the electric-field potential (f satisfy the Laplace equations 

= 0, VV = 
with the following boundary conditions: 

2 Snp p v^l + iV±ri)^ 

r]t = Vl?7 ■ Vl^, z = r]{x, y, t), 
(p = 0, z = r]{x,y,t), 
|V$| ^0, z ^ -oo, 



ip — > —Ez, z oo, (1) 

where a is the surface tension coefficient and p is the density of a medium. 

We are interested in the dynamics of formation of a singular profile for a conducting fluid. 
It is natural to assume that the electric field near the cusp appreciably exceeds the external 
field; i.e., {Vipl ^ E. In this case, the interface evolution is fully determined by the intrinsic 
field, which decreases with distance from the singularity. One can thus use the condition 

iVy^l ^0, z^oo (2) 

instead of the field uniformity condition (1). This agrees with the assumption about the 
universal behavior of a fluid in the formation of a singular surface profile, because it allows 
the fluid motion near the singular point to be analyzed without regard for the particular 
geometry of the problem. The applicability of condition (2) will be discussed below in more 
detail after establishing some regularities for the dynamics of a conducting fluid near the 
singularity. Note that the possibility of secondary Taylor cones nucleating at the already 
formed cones counts in favor of the universal mechanism of formation of conic cusps It is 
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clear from this example that fluid " forgets" the boundary conditions at infinity at the stage 
of collapse. 

Let us consider the most important case of the axially symmetric perturbation of the 
surface. Taking into account that, after substitutions 

r ^ r ATraE""^ , z^zAiraE''^, t ^ t Stt^ aE~^ , 



where r = y/ + y^, the equations of motion become dimensionless and do not contain any 
physical characteristics, one obtains 

^rr + r~^^r + ^zz = 0, z<r](r,t), (3) 
+ r~ Vr- + V'^z = 0, z>r]{r,t), (4) 
+ + + -.(M), (5) 



rit = ^z- Vr^r, z = ri{r, t), (6) 

ip^O, z^r](r,t), (7) 

^l + ^l^O, ^l + ^l^O, r^ + z^^oo, (8) 

$^ = 0, = 0, Vr = 0, r = 0. (9) 
These equations allow the only self-similar substitution 

^x,y,z,t) = ^{f,zy/^ (10) 

cp{x,y,z,t) = (p{f,zy^\ (11) 



v{x,y,t) = fi{f)T 



(12) 
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r = fr^/^, (13) 



z = ~zt''I\ (14) 

where r = tc — t and tc is the collapse time. This substitution occurs due to the fact that 
Eqs. (3)-(9) are invariant about dilatations 

$ — > $c, ip — > ipc, T] —>■ r]c^, 

r — > rc^, 2; — > zc^, t — > te^, 

i.e., it occurs, in fact, from dimensional considerations (c is an arbitrary constant). Note that 
the initial electrohydrodynamic equations with condition (1) do not permit one to introduce 
any self-similar variables. 

Substituting Eqs. (10)-(14) in Eqs. (3)-(9), one finds that the functions l», (p and fj obey 
the following set of partial differential equations: 

^ff + r-^^f + ^zz^O, z< fj{r), (15) 



(ffr + r ^(ff + (fSz^O, z>fj{f), (16) 



3 +2-2 +yrTWVi+^l ^^'^ 



2'qfr - 2r] = - ^f^f: z = r]{r), (18) 



<^ = 0, z^fiif), (19) 



l'? + l>?^0, (^? + (^?^0, + 52^00, (20) 



$- = 0, (ff = 0, r/f= = 0, f = 0. (21) 
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For self-similar solutions (10)-(14), the surface profile forms first at the periphery and then 
extends to the center r = z = (the scale decreases as r^/^). This implies that the formation 
of conic cusps at t = tc is described by those solutions to the set of Eqs. (15)-(21) which 
provide conic asymptotic shape of the surface. In such a situation, the presence of asymptotic 
solutions for which ~ f at f — > oo is the necessary condition for the validity of our 
assumption about the self-similar nature of conic formations. 

Analysis of Eqs. (15)-(21) in the limit R = ^f'^^z^ oo showed that they have an 
asymptotic solution of the form 

<^(f,5) =p-i[2i?(so-s)]^/Vi/2(cos^^), (22) 



<l(f,5) = (23) 

r\{r) = -sof, (24) 
Pi/2(cos^^o) = 0, 

So = -ctg6'o, 

where 6 = arctg(f/5) is the polar distance in spherical coordinates, P1/2 is the Legendre 
polynomial of order i, and s is a constant satisfying inequality < s < sq. This solution 
describes a conic surface with an angle of 2tt — 26q that is equal to approximately 98.6°, i.e., 
to the Taylor cone angle. According to Eq. (23), the fluid motion is spherically symmetric, 
and fluid moves to the sink point R = along the tangent to the surface (24). Since the 
self-similar solution assumes its asymptotic form at r 0, a conic cusp with Taylor angle 
forms at time tc, and Eqs. (22)-(24) are the exact analytic solution of the problem. The 
electric field at the cusp increases as r~^/^, the cusp growth velocity increases as r~^/^, and 
the cusp curvature increases as r~^/^; i.e., these quantities become infinite during a finite 
time. At an appreciable distance from the singularity, the field strength does not change, 
and the velocity of fluid linearly decreases with time and becomes zero at t = tc- The latter 
fact allows one to explain qualitatively the mechanism of transition to the stationary regime 
occurring for liquid-metal sources after the initiation of field ion evaporation from the cusp 
(stationary emitter models were developed in 
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This analysis is only valid on the condition that the asymptotic solutions to the set of 
partial differential Eqs. (15)-(21) have the form of Eqs. (22)-(24). To prove this statement, 
one should construct the asymptotic expansion for the solutions at i? — > oo with leading 
terms given by Eqs. (22)-(24). Let us seek this expansion in the form 

Q'in 

n=0 



^{f,z) = T.^n-Q^[R'^'Pi/2{cos9)] , (25) 



^ir,~z) = Y.bnQ^[R-'], (26) 

n=0 
oo 

f){f)^J2'nr'''''^ (27) 

n=0 

where it is taken into account that the derivative of a harmonic function of any order with 
respect to z is also a harmonic function. The zero-order coefficients are determined by Eqs. 
(22)-(24), 

Oo = p"^[2(so - s)]^/% bo^s, Co = So- 
ft turns out that, to the first order, the surface is conic: 



18so(3 - 2sl) ' 



ai = 0, bi = °r.,2\ ' "^i = 0. 



The correction to Eq. (24) for the surface shape appears in the next order. One finds from 
kinematic boundary condition (18) that 

The coefficients 02 and 62 can be determined from Eq. (19) and, correspondingly, Eq. (17), 
where one should use the linear order of perturbation theory for small deviation of the surface 
from the cone. In turn, the coefficient C3 is found from Eq. (18), where one should take 
into account the quadratic nonlinearity, etc. Thus, the expansion coefficients are uniquely 
determined by the zero-order coefficients, i.e., by the parameter s of the problem, confirming 
the existence of solutions with the desired asymptotic form for the equations of motion. Note 
that, if one sets s = and, hence, $ = 0, then Eqs. (15)-(21) coincide with the Taylor 
equations in the problem of equilibrium configuration of a charged liquid-metal surface. 
However, one fails to construct asymptotic expansion (25)-(27) in this case. 
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FIG. 1: Schematic drawing of a surface of conducting fluid z = rj{r) corresponding to the self- 
similar solutions of the equations of motion. 

Further, the solutions given by expansions (25)- (27) adequately describe the experimental 
data if the condition fj < —s^f is satisfied, i.e., if the fluid surface is positioned above 
the asymptotic cone (see figure P). Otherwise, the surface velocity would be directed in 
opposition to the z axis. Let us check how this condition is fulfilled in the limit of large r. 
It follows from the expansion obtained above for the surface shape that it deviates from the 
conic shape in the direction specified by the sign of the C2 coefficient. In our case, C2 < 
(c2 > for cone angles smaller than 78.5° and larger than 126.9°), so that the amplitude of 
surface perturbation should increase. Indeed, the evolution of the fluid boundary away from 
the singularity is determined by the leading terms of the expansion in small r value, 

//(r, t) = -sqt - |c2|r'^r"^, 

from whence it follows that, when forming a conic cusp, the fluid moves upwards, as is 
expected from physical considerations. 

Let us now consider the surface geometry for small f, where expansion (27) diverges. For 
the function fj to satisfy condition (21), the surface near the cone apex must be "rounded 
off" (figure). Let us estimate the distance |?7(0)| from the cone apex to the fluid surface. 
Multiplying kinematic boundary condition (18) by 27rf/3 and integrating it over f, one 
obtains after simple mathematics 

2V = [ dn^dS, 
Js 

where S stands for the fluid surface z = fj^f), V is the volume of a region bounded from 
above by the conic surface z = — Sor and from below by the S surface, and dn denotes the 
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derivative along the normal to S. The integral on the right-hand side of this expression is 
equal to the fluid velocity flux through the surface. Since the function l» is harmonic, the 
vector-fleld flux Vl* through any closed surface is zero. This fact allows the flux through 

the surface S to be determined using the asymptotic form of velocity potential at i? — > oo. 
Taking into account that the fluid flows into a solid angle 27r(l + cos 6*0) at infinity, one has 
fromEq. (23) 

/ dn^dS ^27rs{l + cos 60), 
Js 

and, hence, V = 7rs(l + cos^^o)- Notice that the volume of a region bounded by the conic 
surface z = —sof and the plane z = —h (a circular right cone of height h) is equal to V at 

h = h{s) = [3sso\l + cos ^0)] ^^'^ ■ 

Clearly, if the volume V is fixed and the conditions r/(f) + sgf < and f/^ (f) < are fulfilled 
for any f, the quantity |f/(0)| cannot exceed the cone height. That is, the inequality 

\m\<h{s), 

connecting the characteristic spatial scale at small it! with the asymptotic parameter s is 
satisfied. Since the most probable value of h{s) corresponds to the maximum allowable value 

So of the s constant, the following estimate is also valid: 

1^7(0)1 </i(so) = So(3 + 3 cos ^o)'/', 

which does not involves the free parameter s. 

Let us return to the question of the applicability of approximation (3)-(9) to the initial 
equations of motion. As was pointed out above, condition (2) can be used instead of Eq. (1) 
only if the external electric field is weaker than the intrinsic cusp field. This implies that the 
inequality (f^ + ifl 1 must be fulfilled. After the transition to the self-similar variables, it 
is recast as 

It is clear that for small r (i.e., immediately before the collapse) this condition is fulfilled 
near the singularity in a natural way. Because one can write |V(^| ~ (r^ -|- z^)~^/^ at small 
T, one has |V</p| ^ 1 in a rather close vicinity of the singularity. In this case, Rq and tq 
values exist for which model (3)-(9) with < r'^ + z'^ < and < r < tq adequately 
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describes the strongly nonlinear stages of electrohydro dynamic instability development for 
the free surface of a conducting fluid in an external electric field. At r"^ + z"^ > Rq, the role 
of nonlinear processes is rather insignificant; the condition for field uniformity at infinity 
(1) should be used together with the corresponding conditions in the limit r — oo, and, 
in particular, with the condition for spatial localization of surface perturbation: rj at 
r oo. In this region, the evolution of fluid surface is described by p erturbation theory for 
small surface slope; this procedure was implemented, e.g., in [l^ 

In conclusion, let us discuss the possibility to form stronger singularities-cuspidal points. 

n 

It is known [12j that the field mainly increases as upon approaching the apex of a thin 
point and, hence, the electrostatic pressure Pe increases as r~^. Since the surface pressure 
changes as Ps ~ r~^, Pe ^ Ps near the singularity, and the capillary effects can be ignored. 
It was shown in that in the absence of surface tension weak root singularities rj ~ r^/^ 
form, for which the curvature is infinite, while the surface itself remains smooth. Therefore, 
when assuming that the cuspidal points may appear, one arrives at a contradiction. This 
gives grounds to assume that conic singularities are precisely those which are the generic 
singular solutions of the electrohydrodynamic equations, so that the behavior of a charged 
liquid-metal surface with cusps is described by self-similar solutions (10)-(21). 

Note also that the results of this work can be extended to dielectric fluids, in which the 
conic cusps with angles depending on the dielectric constant can form in an electric field 
In addition, the approach developed in this work can be applied to the description 
of the evolution of dimples sharpening in a finite time at a liquid helium surface (see, e.g., 
experimental work ^14j). In my preceding work [15] devoted to the construction of exact 
analytic solutions to the equations of motion for liquid helium in the presence of weak 
capillary effects, I proved that cuspidal points r] ~ appear at the surface in planar 

geometry. The question of the singularity type for axial symmetry has not been considered 
so far. 
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